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Abstract. Three-dimensional large eddy simulations of solar surface con- 
vection using realistic model physics are conducted. The thermal structure of 
convective motions into the upper radiative layers of the photosphere, the range 
of convection cell sizes, and the penetration depths of convection are investigated. 
A portion of the solar photosphere and the upper layers of the convection zone, 
a region extending 60 x 60 Mm horizontally from Mm down to 20 Mm be- 
low the visible surface, is considered. We start from a realistic initial model of 
the Sun with an equation of state and opacities of stellar matter. The equa- 
tions of fully compressible radiation hydrodynamics with dynamical viscosity 
and gravity are solved. We use: 1) a high order conservative TVD scheme for 
the hydrodynamics, 2) the diffusion approximation for the radiative transfer, 
3) dynamical viscosity from subgrid scale modeling. The simulations are con- 
ducted on a uniform horizontal grid of 600 x 600, with 168 nonuniformly spaced 
vertical grid points, on 144 processors with distributed memory multiprocessors 
on supercomputer MVS-15000BM in the Computational Centre of the Russian 
Academy of Sciences. 



Introduction 



Convection near the solar surface has a strongly non-local and dynamical char- 
acter. Hence, numerical simulations provide useful information on the spatial 
structures resulting from convection and help in constructing consistent models 
of the physical processes underlying the observed solar phenomena. We conduct 
an investigation of the temporal evolution and growth of convective modes on 
scales of mesogranulation and supergranulatio n in a thre e -dime nsional compu- 
tational box. In previous work by the author Ustvugov ( 20061 )] it was shown 



that collective motion of small convective cells of granulation expels weak mag- 
netic field on the edges of cells at mesogranular scales. The average size of such 
cells is 15-20 Mm and th e lifetime of order 8 -10 so lar hours. Simulation of solar 
photosphere convection Stein Nordlund (2006)] in a computational domain 



of size 48 Mm in the horizontal plane and 20 Mm in depth showed that the sizes 
of convective cells increase with depth. The purpose of this work is to investi- 
gate the development and scales of convection in a region of size 60 Mm in the 
horizontal plane and 20 Mm in depth. 



1 



2 



Ustyugov 



Numerical method 



The distribution of the main t hermodynamic variable s with radius was taken 
from the Standard Solar Model [Christensen-Dalsgaard ( 20031 )] with parameters 
(X,Z,a) = (0.7385,0.0181,2.02), where X and Y are hydrogen and helium 
abundances by mass, and a is the ratio of mixing length to pressure scale height 
in the convection regi on. We used the OPAL o pacity tables and the equation of 
state for solar matter [iglesias Sz Rogers! ( 19961 )]. 

The fully compressible nonideal hydrodynamics equations were solved: 



dp 
dt 



+ V • pv = 



(1) 



dpv 

~~dT 



+ V 



pvv + PI = pg + V • f 



(2) 



dE 
~dt 



+ V ■ [v(E + P)] = V ■ (v ■ f) + p(g ■ v) + Q rad 



(3) 



E = e + pv 2 /2 is the total energy, Q ra d is the energy transferred by radiation and 
f is the viscous stress tensor. The influence of small scales on large scale flows 
was evaluated in terms of the viscous stress tensor and the rate of dissipa t ion w as 



defined from the buoyancy and shear production terms [Canuto et.al.1 (|l994l )]. 
The evolution of all variables in time was found usin g an expli c it T VP [Total 
Variation Diminishing] conservative difference scheme Yee et.al. (jl990l )]: 
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The flux along each direction was defined by the local-characteristic method as 



i+l/2,j,k 



Ei,j,k Pi 



i+l,j,k 



4+1/2^1+1/2 



(6) 



R-i+i/2 is the matrix whose columns are right eigenvectors of dF/dU evaluated 
as a generalized Roe average of and for real gases. VFi+1/2 is the 

matrix of numerical dissipation. The term Sij t k represents the effect of gravita- 
tional fo rces and radiation. T ime step integration is by third-order Runge-Kutta 
method Shu Osherl ( 1988)]. This scheme is second-order in space and time. 
Central differences were used for the viscous term, and the diffusion approxima- 
tion was applied for the radiative term in the energy equation. 

We used a uniform grid in the x and y directions and a nonuniform grid 
in the vertical (z) direction. Periodic boundary conditions were used in the 
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horizontal directions and the top and bottom boundary conditions were chosen 
to be 

v z ,k = ~v z ,k-i,v x ,k = v x ,k-i,v y ,k = v y ,k-i,dp/dz = pg z ,p = p(p),e = const 

that is, reflection for the z component of velocity, outflow for the x and y com- 
ponents. Pressure and density were derived from the solution of the hydrostatic 
equation, using the equation of state with a constant value of the internal energy. 

3. Results 

The numerical simulation of solar photosphere convection was conducted during 
24 solar hours. On the upper boundary of the region, which corresponds to op- 
tical depth r = 1 near the solar surface in the horizontal plane, the development 
of small-scale convection structure is clearly seen, with cells of average size 1.5 
Mm and lifetime of order 1-2 Min (Figure 1). At depth 1 Mm from the solar 
surface we observe the appearance of signs of large-scale convective cells at the 
scale of mesogranulation (Figure 2). With increasing depth the size of big con- 
vective cells grows and their number decreases. At depth 2-5 Mm the average 
size of a cell is 10-15 Mm, whereas it is 15-20 Mm at depths greater than 10 Mm 
(Figure 3-6) . Vertical motion on the edges of convective cells has the maximum 
values of velocity and become supersonic at cell vertices at depths 3-4 Mm. In 
the deepest part of the domain the flow of material becomes slower and more 
laminar. 

By looking at the convective flow in a vertical plane we can see three re- 
gions, distinctly different in physical character (Figure 7-8). In the turbulent 
zone at depths 0-2 Mm, small-scale chaotic motion of granules due to the action 
of strong compressible effects leads to formation of downdrafts of cold blobs of 
material. There are places with strong vorticity motion where material moves 
from different directions to one point with the formation of more powerful down- 
drafts. In the transition zone at depths of 2-5 Mm we find large-scale flow, on the 
scales of mesogranulation. Here the downdrafts achieve a maximum velocity, on 
average, of about 4 km/sec with Mach number 1.5. Material moves downwards 
in the form of narrow jets and upwards in wider regions with an average of size 
10-15 Mm. In the third part of the region, at greater depths, we obesrve some 
points where material reaches the bottom boundary with subsonic velocity. The 
distances between these points gives us the size of convective cells on the scale 
of supergranulation, about 20 Mm. 

We averaged the velocity components over 4 solar hours for one horizontal 
plane near the solar surface. We then introduced many Lagrange particles or 
corks for each cell, uniformly distributed across the plane at some moment of 
time. After 4 solar hours of simulation we defined the places with the maximum 
number of corks. From Figure 9 it is clearly seen that large-scale divergent flows 
from the centers of big convective cells expel corks to their edges. Corks are 
located in places with converging flows (Figure 9) where we have positive values 
of the two-dimensional divergence of velocity and at points with strong vorticity 
motion where we find large negative values of curl velocity. As seen from Figure 
10 there are convective cells with the size of supergranulation - about 20 Mm in 
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diameter. Matter move away from the center of these cells with average velocity 
1-1.5 km/sec. 
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Figure 2. Contours of the vertical component of velocity at depth 1 Mm. 
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Figure 3. Contours of the vertical component of velocity at depth 2 Mm. 
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Figure 6. Contours of the vertical component of velocity at depth 10 Mm. 
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Figure 8. Contours of vertical component of velocity field. 
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Figure 10. Null contour of two-dimesional divergence of velocity and bound- 
aries of supergranulation cells. 



